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Using squeezed states it is possible to surpass 
the standard quantum limit of measurement un- 
certainty by reducing the measurement uncer- 
tainty of one property at the expense of another 
complementary property [Ij. Squeezed states 
were first demonstrated in optical fields [2j and 
later with ensembles of pseudo spin-1/2 atoms 
using non-linear atom-light interactions [3j. Re- 
cently, collisional interactions in ultracold atomic 
gases have been used to generate a large degree 
of quadrature spin squeezing in two-component 
Bose condensates |4l[5]. For pseudo spin-1/2 sys- 
tems, the complementary properties are the dif- 
ferent components of the total spin vector (S), 
which fully characterize the state on an SU(2) 
Bloch sphere. Here, we measure squeezing in a 
spin-1 Bose condensate, an SU(3) system, which 
requires measurement of the rank- 2 nematic or 
quadrupole tensor (Qij) as well to fully charac- 
terize the state. Following a quench through a 
nematic to ferromagnetic quantum phase tran- 
sition, squeezing is observed in the variance of 
the quadratures up to -8.3^0 ? dB (-10.3^0.9 dB 
corrected for detection noise) below the standard 
quantum limit. This spin-nematic squeezing is 
observed for negligible occupation of the squeezed 
modes and is analogous to optical two-mode vac- 
uum squeezing. This work has potential applica- 
tions to continuous variable quantum information 
and quantum-enhanced magnetometry. 

The study of many-body quantum entangled states in- 
cluding atomic spin squeezed states is an active research 
frontier. In addition to being intrinsically fascinating, 
such states have applications in precision measurements 
[6^ , quantum information and fundamental tests of quan- 
tum mechanics [7 . Atomic squeezed states were first 
considered for ensembles of two-level (pseudo spin-1/2) 
atoms. For spin-1/2 particles, coherent states of the sys- 
tem are uniquely specified by the components of the total 
spin vector (S), typically illustrated on a SU(2) Bloch 
sphere. For particles with higher spin, additional de- 
grees of freedom beyond the spin vector are required to 
fully specify the state. For spin-1 particles, a natural ba- 
sis to describe the wavefunction is the SU(3) Cartesian 
dipole-quadrupole basis, consisting of the three compo- 
nents of the spin vector, 5^, and the moments of the rank- 
2 quadrupole or nematic tensor, Qij {{i^j} G {x^y^z}). 
In matrix form, the nematic moments can be written 
Qij — SiSj + SjSi — -^Sij [8 . 

Spin-1 atomic Bose-Einstein condensates [9HT3] pro- 
vide an exceptionally clean experimental platform to in- 
vestigate the quantum dynamics of many-body spin sys- 



tems. They feature controllable quantum phase tran- 
sitions, well-understood underlying microscopic models, 
and flexible defect-free geometries. Importantly, it is pos- 
sible to initialize non-equilibrium or excited states of the 
system and to directly measure both the spin vector and 
the nematic tensor using standard atomic state manipu- 
lation tools. Law, et al.^ demonstrated that the spinor in- 
teraction can be written as total spin angular momentum, 
A^^ where S'^ = Sl^Sl^Sl [14 . It is natural to describe 
this many-body system in the second quantized formal- 
ism in terms of the mode operators of three Zeeman states 
(e.g. Sx = -^{a\ao + aj^-i + + ^-i^o))- Here 
are particle annihilation operators, and we use the sin- 
gle mode approximation where the modes share the same 
same spatial wavefunction 0(r) determined by the spin- 
independent part of the Hamiltonian. The spinor Hamil- 
tonian describing the collisionally-induced spin dynamics 
of the condensate and the effects of an applied magnetic 
field B (by convention along the z-axis) can be written 
(Supplementary Information): 

H = XS^^ ^qQ,, (1) 

Here A and q (x characterize the inter-spin and Zee- 
man energies, respectively. At high fields, the system fa- 
vors nematic ordering of the spins, an eigenstate of Qzz- 
This is a state with (S) = with broken rotational sym- 
metry given by the anisotropy of spin fiuctuations, e.g. 
{SD = (Sy) 7^ (5^) and whose alignment is specified 
by a time-reversal invariant 'director', in this case the 
z-axis. In the Fock basis, |A^i, A^o, ^-i), this is just the 
state with all N atoms in the mj = state: |0, 0). At 
low fields, the sign of A determines whether the interac- 
tions favor a ferromagnetic (A < 0, i.e. ^^Rb used in our 
work) or anti-ferromagnetic (A > 0, i.e. ^^Na) ground 
state with the ground states being maximum and mini- 
mum total spin respectively. At intermediate fields, the 
system undergoes a quantum phase transition between 
orders with a quantum critical point at q = 2\c\ for the 
ferromagnetic case, where c = 2NX [15 . 

In spin-1/2 systems, non-linear Hamiltonians such 
as S'^ produce squeezing of the spin variables satisfy- 
ing Heisenberg uncertainty relations, e.g. ASx^Sy > 
While criteria for squeezing and entanglement 
are now well-established for spin-1/2 particles within a 
SU(2) framework, there was considerable early discussion 
about different squeezing definitions [16, 17 . There has 
been much less work for higher spin particles with corre- 
spondingly higher symmetries. Squeezing in spin-1 sys- 
tems has been studied from the perspective of two-mode 
squeezing |I8H2T] . in terms of the Gell-Mann (quark) 
framework of the SU(3) algebra [22^, and in terms of 
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FIG. 1: Illustration of the experimental sequence using semi-classical simulation and quasi-probability distributions, a. The initial state 
is a condensate with the atoms prepared in the ruf = state. A N = 30 atom distribution is used to emphasize features, b. After 25 ms 
of evolution, spin-nematic squeezing develops along the separatrix (green line) in the two upper spheres, c. A microwave pulse rotates 
the quadrature phase. For comparison the state from the previous plot is shown in red in the upper two spheres, d. A 7r/2 RF pulse 
rotates the transverse magnetization Sx into Sz- For comparison the state from the previous plot is shown in red in the lower two spheres, 
e. After the trap is turned off, a Stern-Gerlach field is applied during the TOF expansion and the clouds of atoms are counted using 
fiuorescence imaging. Shown are points for 100 runs of a squeezed quadrature (green) and an unsqueezed quadrature (orange) from the 
experimental measurements. 



spin-nematic measurables [8]. Although in each study 
appropriate phase spaces for squeezing were identified, a 
comprehensive picture remains to be developed. 

The source of squeezing in a spin-1 condensate is the 
non-linear collisional spin interaction Tig = AS^, which 
reduces to 5^ + Sy for states with (Sz) =0 of in- 

terest here. This Hamiltonian contains four- wave mixing 
terms, 'Hfwm = 2A(aQ^aia_i + aJaLi^o)? which gener- 
ate squeezing that can be described using a two-mode 
formalism where the a±i modes are the signal and idler 
(Supplementary Information). Here we prefer an analysis 
using the commutators of the SU(3) group describing the 
spin-1 system as it provides more insight. Experimental 
investigations of spin-1 condensates have been predomi- 
nantly in the mean-field limit, however first explorations 
beyond the mean field have been reported [23H26] . 

From the generalized uncertainty relation AAAB > 
^])|, only operator pairs with non-zero expecta- 
tion values for their commutation relations can exhibit 
squeezing. For condensates with the atoms in the m/ = 
state, only two of the SU(3) commutators have non- 
zero expectation values : (0, A^, 0| [^S^^, Qyz] |0, A^, 0) = 
-2iN and (0, AT, 0| [S'^, Q^^] |0, A", 0) = 2iN, leading to 



the relevant uncertainty relations ASxAQy^ > A" and 
ASyAQxz ^ A'. For each of these, the uncertainty re- 
lationship is between a spin operator and a quadrupole 
nematic operator; these operators and their commutators 
define two subspaces. From these relations, two squeez- 
ing parameters are defined in terms of quadratures of the 
operators: 

e^^y^ = {(A(cos^5,(,) + sin^Q,,(,,)))2)/Ar. (2) 

with 6 as the quadrature angle [3 [20] . Squeezing within 
a given SU(2) subspace is indicated by the variance of the 
quadrature operator being less than the standard quan- 
tum limit (SQL) of A^ for some value of 6. 

The experimental sequence is illustrated in Fig. 1 with 
the help of the spin SU(2) subspace {Sx^ Sy^ Sz} and both 
of the subspaces that exhibit squeezing, {Sx^Qyz^Qzz} 
and {Sy^Qxz^Qzz} (Supplementary Information). The 
initial state of the condensate is shown in Fig. la. It 
has no spin moment, (S) = 0, but non-zero diagonal 
quadrupole elements (Qu) ^ 0, and uncorrelated equal 
uncertainties in Sx, Sy^Qxz^Qyz- Subsequent evolution 
is governed by Eq. 1 in the regime q < 2\c\. Out of 
equilibrium dynamics of spin-1 condensates generally ex- 
hibit oscillatory behavior in the spin components [27l [28] , 
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FIG. 2: Comparison of measured quadrature fluctuations with a fully quantum calculation, a. Measurement of the quadrature variances 
for different evolution times and quadrature phase shifts. The phase shift is calculated from the microwave detuning normalized to the 
on-resonance Rabi rate (Supplementary Information). Estimated errors are approximately the size of the marker for both phase and 
variance. Open markers are statistics of the raw data, fllled markers have been corrected for PSN. b. Detailed view of the maximum 
squeezing for different evolution times. The phase of data is shifted to match the simulation, c. Time evolution of the populations and 
the maximum and minimum quadrature variances. The squeezing measurment is performed prior to signiflcant population evolution. 



except near a separatrix (green contour) in phase space 
where the period diverges [29]. For our case, the ini- 
tial state lies on the separatrix, and hence, evolution is 
wholly dictated by quantum fluctuations corresponding 
to the initial uncertainties. These uncertainties evolve 
exponentially, generating anti-squeezing for a quadrature 
of {Sx^Qyz} aligned along a branch of the separatrix 
and squeezing for the orthogonal quadrature, as shown in 
Fig lb. Measuring the squeezing requires state tomogra- 
phy involving two SU(3) rotations. The first is a rotation 
about Qzz that rotates the quadrature squeezing to align 
to Sx' Conceptually, Qzz is the generator of the rota- 
tion in the squeezed subspaces and the rotation of the 
quadrature angle is given by the operator ex.p{i A6 Q^z)- 
The second is a 7r/2 rotation about the Sy axis (in the 
lab frame) that rotates the fluctuations in Sx into the 
measurement basis, Sz = Ni — N-i (Fig. Id). These 
manipulations are equivalent to the homodyne technique 
used in quantum optics to measure two-mode squeezing 
where the do mode is the local oscillator (Supplemental 
Information). 

Identical squeezing also occurs in the degenerate 
{Sy^Qxz-tQzz} subspace, which leads to an important, 
but subtle point. In the lab frame, the system does not 
confine itself to either of the SU(2) subspaces but rather 
undergoes rapid Larmor precession described by rota- 
tions about Sz in two other SU(2) subspaces {Sx-, Sy^ Sz} 



and {Qyz^Qxz^ Sz} (not shown). However, because the 
squeezing is identical in both subspaces, and the Larmor 
precession of the spin vector and quadrapole are synchro- 
nized, it is not necessary to track the precession in order 
to measure squeezing. This has important experimental 
consequences in that it is not necessary to maintain syn- 
chronization with the Larmor rotation in order to perfom 
quantum state tomography. 

For the experiment, we prepare a condensate of = 
45, 000 ^^Rb atoms in the |/ = 1, mj = 0) hyperfine state 
in a high magnetic field (2 G). The condensate is tightly 
confined in an optical dipole trap with trap frequencies of 
250 Hz. For these parameters the Thomas- Fermi radii of 
the condensate are 3.8 jam while the spin healing length 
is 11 jam indicating that the condensate is well described 
by the single mode approximation. To initiate evolution, 
the condensate is quenched below the quantum critical 
point by lowering the magnetic field to a value 210 mG 
and then allowed to freely evolve for a set time. During 
this evolution is when the squeezing is developed. Note 
during this time there is essentially no population trans- 
fer (< 1%) from the rrif = state (Fig.|2j3 top), hence this 
corresponds to squeezed vacuum of the m/ = ±1 modes. 
The quadrature rotation is accomplished by shifting the 
phase of the m/ = amplitude using a microwave de- 
tuned from the clock transition. This is immediately fol- 
lowed by an RF rotation to rotate Sx into Sz- The trap 
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FIG. 3: Reconstructions of the phase space for different evolution times: a 15 ms, b 30 ms, c 45 ms, and d 65 ms. The last two are at 
half the scale of the first. The black trace in each is the calculated 1/Ve uncertainty ellipse from the simulation. 



is then turned off and a Stern-Gerlach field is applied to 
separate the rrif components during 22 ms time-of- flight 
expansion. The atoms are probed for 400 /is with three 
pairs of orthogonal laser beams, and the resulting fluo- 
rescence signal is collected by a CCD camera. Sample 
Sz measurements are shown in Fig. le for two different 
quadrature angles and show qualitatively both enhanced 
and suppressed fluctuations. 

The experiment cycle is repeated 100 times for each 
quadrature angle to collect sufficient statistics to deter- 
mine the fluctuations The measurement results are 
shown in Fig. [2] for different evolution times. The fluctua- 
tions clearly exhibit noise reduced below the SQL (^^ < 
dB) for certain quadrature phases and show increased 
noise 7r/2 radians away. The measurements are com- 
pared to a fully quantum theoretical calculation (Supple- 
mental Information) performed by numerically integrat- 
ing Eq. ([T]) with a spinor dynamical rate c/h = —8 Hz 
{h is Planck's constant) and q{t) determined experimen- 
tally using microwave spectroscopy. The spinor dynam- 
ical rate is chosen to match the degree of anti-squeezing 
observed and is also a good fit to the long time evolution 
of the populations (Fig.[2]3 top) as well as estimates from 
the trap frequencies and N. 

As evolution time is increased, the maximum variance 
shows anti-squeezing that increases in good agreement 
with the simulation (Fig.|2]3. black lines) and is approx- 
imately exponential with a time constant of (2lc|/fi)~^ 
given by the two-mode squeezing model (Fig.|2]3. thin 
gray lines). The minimum variance shows squeezing that 
limits asymptotically after ~ 30 ms of evolution due to 



detection noise from a combination of light scattered by 
the apparatus and the photo-electron shot noise (PSN). 
The PSN limit is indicated by the gray lines in Fig.[2j 
The maximum observed squeezing is — 8.3^'^Q y dB, which 
is the highest degree of quadrature squeezing observed 
in any atomic system. When corrected for the PSN, 
it is possible to infer a 'corrected' squeezing parameter 
of — 10.3lo!9 ^^^^ would be obtained with detection 
improvements. The phase of maximum squeezing also 
evolves in time, converging to the phase of the separatrix 
given by cos^s = —q/c — 1 with a small (~ 150 mrad) 
discrepancy between the measured phase of maximum 
squeezing and the theoretical prediction. After this ini- 
tial time of approximately exponential squeezing, the 
minimum squeezing parameter reaches a turning point 
due to pump depletion. 

We also reconstruct the phase space distribution of the 
squeezing for each time (Fig.[3| using an inverse Radon 
transform. (Supplementary Information). The theoret- 
ical squeezing ellipse is shown for comparison. Just as 
the theoretical prediction shows the quadrature evolves 
both in aspect ratio and phase rotation, however, quan- 
titative agreement is limited by the finite number of data 
sets. Furthermore the minimum width of parts c and d 
is limited by the PSN. 

In conclusion, we have observed spin-nematic squeezed 
states created from the free-evolution of a spin-1 conden- 
sate after being quenched through a nematic-ferromag- 
netic quantum phase transition. This work demonstrates 
quantum dynamics of a mesoscopic SU(3) spin-nematic 
many-body system and paves the way for further investi- 
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gations beyond the mean- field limit. It shows that simple 
evolution of a spin-1 condensate is a robust way to gen- 
erate spin-nematic squeezed states with large amounts 
of squeezing. These states, together with the tools we 
have demonstrated for their characterization, can be used 
for quantum metrology of magnetic fields and for atomic 
clocks. These experiments demonstrate new paths to 
explore the fascinating intersections of correlations, en- 
tanglement and quantum phase transitions in an exotic 
quantum magnetic system. 

Very recently, two papers related to ours have appeared 
[30, 31 . These experiments share a similar spin mixing 
mechanism to generate quantum correlations, though in 
their case it is in the / = 2 hyperfine manifold. For 
both of these other experiments, the starting point is 
the observation of strong non-classical fiuctuations of 
Sz = J^i — J^-i due to spin mixing, which we previously 
reported in Ref. [26 . In Ref. [31 , atomic homodyne de- 
tection of the quadrature fiuctuations using a technique 



that is very similar to ours was demonstrated. In their 
case, they were not able to show that the fiuctuations 
were smaller than the classical limit, only that the fiuc- 
tuations were dependent on the local oscillator phase. 
Only by subtracting off detection noise do they claim en- 
tanglement of the two modes at 4 ± 17 atom level. In 
Ref. [30 , the number-squeezed mj = ±1 states are cou- 
pled using a sequence of microwave transitions to realize 
a two- level pseudo spin-^ system. On the corresponding 
Bloch sphere, the number squeezed states are completely 
phase uncorrelated. By measuring the resulting variance 
and 4*^ moment of Sz, they were able to show that the 
angle could be determined with an uncertainty -1.6 dB 
below the classical limit. These systems are in the regime 
of twin-atom states similar where the mj = ±1 states 
are macroscopically populated. In this regime there is no 
longer a simple spin-nematic squeezed vacuum interpre- 
tation of the squeezing, but it is still possible to measure 
quantum correlations as these experiments demonstrate. 
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